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Coherence of spins of electrons confined in III-V quantum dots is strongly affected by their hy- 
perfine interaction with the nuclei. In this paper an introduction to this subject is presented. 
Some theoretical approaches to the problem will be outlined. Most attention will be given to the 
Quasi-Static Bath Approximation, to the cluster expansion theories of dephasing due to the nuclear 
dynamics induced by the dipolar interactions (spectral diffusion), and to the effective Hamiltonian 
based theory of dephasing due to hyperfine-mediated interactions. The connections between the 
theoretical results and various experiments will be emphasized. 
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1. Introduction 

Coherently controlled electron spins embedded in a 
semiconductor environment have been a subject of in- 
tense research during the last ten years. The main mo- 
tivation for the experimental and theoretical work was 
the prospect of their use as qubits [T] , the basic building 
blocks of quantum computers. Recently, coherent control 
of a single spin was also shown to allow for very sensitive 
detection of magnetic fields with nanoscale resolution [2] . 

During the last few years a large progress was made 
in coherent control of spins in III-V materials such as 
GaAs [5HTU] and InGaAs [TTJ [T3] . In these semiconduc- 
tors all the nuclei are spinful, and the hyperfine (hf ) in- 
teraction with the nuclear spins is the main factor af- 
fecting the electron spin coherence time. Due to a large 
mismatch of Zeeman energies between the electron and 
the nuclei, the latter affect very little the longitudinal 
electron spin relaxation, which is thus due to phonon- 
induced transitions between the electronic states having 
a mixed spin character because of spin-orbit interaction, 
see e.g. [HESE]. Let us note that these interactions lead 
to spin relaxation (and spin dephasing |14j ) on timescales 
longer than miliseconds [4 , which are much longer than 
the timescales considered here. The hf interaction with 
the nuclei is however the dominant source of the trans- 
verse spin decay, i.e. the electron spin dephasing. Since 
the preservation of a specific phase relation (coherence) 
between the two states of the spin-1/2 is the necessary 
condition for the spin to be considered a qubit, the hf- 
induced decoherence needs to be understood. 

The aim of this paper is to provide an introduction to 
the problem of hf-induced electron spin decoherence. In 
contrast to many other problems which involve decay of 
a population of a quantum state or randomization of a 
relative phase between two such states (e.g. spontaneous 
recombination and dephasing of optical transitions), it is 
impossible use the Bloch equations [HI [16] to describe 
the dynamics of the electron spin interacting with the 
nuclei. In order to see this one should remember the 



conditions under which these equations are derived: the 
coupling between the system of interest (the electron spin 
here) and the bath should be weak (i.e. the interaction 
should be small compared to the self-Hamiltonians of the 
spin and the bath) , and the autocorrelation time of the 
bath should be much shorter than the timescale of the 
system's dynamics. Only then one obtains that the off- 
diagonal elements of the density matrix of the system 
decay according to p a b(t) ~ exp(— t/T-z). This is not the 
case for the electron spin interacting with the nuclei: the 
nuclear dynamics is very slow, and the coupling of the 
spin to the nuclei can be considered weak only in very 
high magnetic fields, B>1 T, and only at these fields 
(and for a specially prepared state of the nuclei) the ex- 
ponential decay of the transverse spin has been predicted 

USE!]. 

We will attempt to give an overview of theoretical ap- 
proaches to this problem, with special focus on theories 
desribing the pure dephasing situation, in which the lon- 
gitudinal relaxation of the spin is neglected, and only the 
dephasing (the decay of the spin component transverse 
to an applied magnetic field) is considered. These are the 
theories of decoherence due to the dipolarly-induced dy- 
namics of the nuclei (the so-called spectral diffusion) [2"2T - 
[2"5] , and the theories focusing on hf-induced dynamics de- 
scribed by an effective Hamiltonian (containing so-called 
hyperfine-mediated interactions) [13 QH [20l EU [26 ] . The 
predictions of these theories for the spin echo decay in a 
GaAs QD have been recently confirmed experimentally 

The only analytical approach not manifestly of the 
pure dephasing kind which will be discussed here is 
the Quasi-Static Bath Approximation (QSBA) [29H3i"j . 
which is applicable at the timescale at which the nuclear 
dynamics is absent or can be treated in a very simplified 
manner (i.e. by replacing all the nuclear spins by a col- 
lective classical vector). There are other approaches to 
the decoherence problem which do not use the effective 
Hamiltonian, and with which the long-time decoherence 
was calculated [STJ I52H54] . The relation between these 
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theories and the effective Hamiltonian based approach is 
a subject of current research. However, since the spin 
echo signal has not been calculated using these methods, 
and since we want to emphasize here the connections be- 
tween the theory and current experiments, we will not 
discuss these approaches further. The exact numerical 
studies of spin decoherence will also be mentioned rather 
briefly, and an interested reader is referred to a review 
[35]. 

Let us mention what else is not covered in this paper. 
Most importantly, we do not explain how the spin initial- 
ization, manipulation, and readout are done experimen- 
tally. The III-V based spin qubits are controlled either 
electrically (by changing the voltages on the gates defin- 
ing a QD in a two-dimensional electron gas), or optically 
(by applying short, coherent, and properly shaped pulses 
of light resonant with optical transitions in self-assembled 
QDs). The physics of each of the control schemes is a 
huge topic in itself, and the interested reader is referred 
to review papers (e.g. [I] for electrically controlled QDs, 
and [351 [37] for the optically controlled ones) . A related 
omission is the lack of discussion of intricate physics of 
the singlet-triplet qubit realized in gated QDs (in this 
qubit the two states of interest are the singlet and un- 
polarized triplet states of two electrons residing in two 
coupled QDs). For simplicity, we focus on the case of a 
single spin interacting with a nuclear bath, but we note 
that all of the results on decoherence are, after suitable 
modifications, applicable to the singlet-triplet qubit (see 
e.g. [35] for the case of Rabi oscillations and [5J for the 
case of spin echo decay) . Finally, let us mention that re- 
cently spins of heavy holes confined in QDs have started 
to attract attention. The interesting features of hyper- 
fine coupling between hole spins and nuclear spins are 
discussed in [351110]. 

We focus on III-V based QDs, but the theoretical ap- 
proaches discussed here are also applicable (with some 
modifications) to other systems in which a localized elec- 
tronic spin interacts with the nuclei. Most importantly, 
the theory of decoherence due to the nuclear dynamics 
induced by the dipolar interactions (discussed in Section 
5) was applied to the case of an electron bound to a phos- 
phorus impurity in Si E3J EH EH ETJ [42] and it shows an 
excellent agreement with the experimental results on spin 
echo in this system j33J [H] . Let us note here that in the 
case of Si one can get rid of nuclear spins by means of 
isotopic purification (removing 29 Si and leaving only the 
spinless 28 Si). The coherence is then limited by inter- 
action with remote electron spins of other donors [42], 
or with the spins of the dangling bonds on Si/SiC>2 in- 
terface [45] . Such purification is impossible in GaAs or 
InAs, since there are no spinless isotopes of Ga, As, and 
In. 

The paper is organized in the following way. In Section 
2 we introduce the Hamiltonian of the system (electron + 
nuclear bath), discuss the characterstic energy- and time- 



scales, and sketch the derivation of the effective Hamilto- 
nian which is applicable at not-very-low magnetic fields 
(we believe it is applicable down to ~ 10 mT in GaAs 
dots, at least at timescales of current experimental inter- 
est). In Section 3 we briefly discuss the experiments in 
which the spin decoherence time is measured. Motivated 
by the slowness of the bath dynamics (mentioned in Sec- 
tion 2), in Section 4 we present theoretical results for 
electron spin dynamics obtained within the QSBA. The 
spin dephasing due to dipolarly-driven nuclear dynamics 
is then discussed in Section 5. Finally, the theory of spin 
decoherence due to hf-mediated interactions (which are 
the dominant sources of decoherence at small magnetic 
fields) is outlined in Section 6. 



2. Basic facts about the nuclear bath 

The Hamiltonian of the electron spin interacting with 
the nuclear bath is given by 



H — Hz + -ffdip + H) 



hf 



(1) 



with the terms corresponding to Zeeman, dipolar, and hf 
interactions, respectively. The Zeeman term is written 
for the magnetic field B in the z direction as 



H 7 , 



UJ. 
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I z 



(2) 



with the electron spin splitting f2 and nuclear spin split- 
tings uj a = —~f a B for distinct nuclear species (i.e. distinct 
elements or isotopes having gyromagnetic ratios 7 Q ), and 
with i being the site index (the nuclear species index a[i] 
is assigned randomly to possible sites while maintaining 
a given ratio of concentration of species) . In the dipolar 
term we employ the secular approximation |46j . i.e. we 
keep only the interactions which conserve the Zeeman 
energy: 



(4 



(3) 



where the summation is over the nuclei i and j of the 
same species, and the couplings given by 



1 1-3 cos 2 0« 

-r 



(4) 



where is the distance between the two nuclei and 
is the angle of relative to the B field direction. For 
nearest neighbors we have 6y~0.1 neV in GaAs (which 
corresponds to the timescale of h/bij~ 10 ms). Together 
with the fact that the nuclear Zeeman energies of Ga 
and As nuclei at B = 1 T are of the order of tens of neV 
(corresponding to less than a milikelvin) , this leads to the 
conclusion that at typical experimental temperatures T 
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we have {Hz + -ffdi p )/fc,BT<C 1. The nuclear bath is thus 
described by a high-temperature density matrix /5b ~ 1- 
The hyperfine Hamiltonian is i?hf = X)i AiS-Ji, which 
we will write as H^f = VqS z + Vg, where we have defined 
the Overhauser operator 

Vb=£>#, (5) 

i 

and the electron-nuclei flip-flop operator 

V^=\Y. A ^ + J- +S-J+) . (6) 

i 

These two are playing very different roles in the process 
of electron decoherence at finite magnetic field. 

In the above equations the hf couplings are Ai = 
where fa = \^{ri)\ 2 , the squared modulus of 
the electron envelope function at the i-th nuclear site 
(with normalization to the primitive unit cell volume: 
J v |^(r)| 2 dr = uq). The hf energies for a nuclear species 
a are A a = |^o^ 2 7s7Jq| u q| 2 : where hq is the vacuum 
permeability, 7s and 7j q are the electron and nuclear 
spin gyromagnetic factors, respectively, and u a is the am- 
plitude of the periodic part of the Bloch function at the 
position of the nucleus of a species (the normalization is 
/ |w(r)| 2 dr = l). For Ga, As, and In atoms A a & 30-50 
fieV (see e.g. [20] or [27] and references therein). The 
largest hf coupling is A max ^A/N, with N being the ef- 
fective number of nuclei in the dot (the number of the 
nuclei which are appreciably coupled to the electron). 
We use a common definition of N = J ^ /i fi/Yli Sii which 
leads to a relation J2iea = n aA 2 t /N, where n a is the 
number (per unit cell) of nuclei of a species. 

It is crucial to notice that due to a large difference 
in electronic and nuclear magnetons, in a finite B field 
we have fl ~ 10 3 uj a . Because of this Zeeman energy 
mismatch a flip-flop between the electron and a nucleus 
is practically prohibited. It can only occur when the 
magnetic field is so low that ft becomes comparable to the 
dipolarly-broadened nuclear lincwidth. This linewidth in 
III-V materials is of the order of a few kHz [35], which 
corresponds to magnetic fields of ~ 0.1 mT. At much 
larger fields the direct electron-nuclear spin flip-flop is 
forbidden, and one can take into account the effect of Vg 
term perturbatively. 

In the second order of perturbation theory we have 
processes in which the electron spin flip-flops with the 
i-th nuclear spin (leading to a virtual state with energy 
differing by ps from the energy of the intial state), 
and then flip-flops back with the j-th nuclear spin. This 
leads to appearance of an effective hyperfine-mediated in- 
teraction between the nuclei, which can be viewed as an 
analogue of the well-known RKKY interaction, only de- 
rived using a localized electronic wavefunction instead of 
Bloch waves. Formally, this interaction is derived by per- 
forming an approximate canonical transformation. The 



effective Hamiltonian is H = e s He s with an unitary op- 
erator e~ s chosen in such a way that Vg is removed from 
the Hamiltonian (for more details see [15 ] 120 ] 125 ] 149 ] ). 
It should be noted that the transformation of states, 
|0) = e~ s \(f)), which should accompany the transforma- 
tion of the Hamiltonian, is not performed in the theories 
of decoherence which are employing the effective Hamil- 
tonian. This approximation is the price which has to be 
paid for obtaining a theory with a convenient structure. 

In the lowest order in Vg we have 5~ VgS z /Sl, and the 
resulting effective Hamiltonian contains the terms 

f" = !(•<?- «>=)+ 

i i 

^"L^tij. P) 

In this Equation the last term is the most important: it 
is the hf-mediated interaction. As expected, its strength 
is decreasing with increasing magnetic field. The crucial 
feature of this interaction is that it is long-ranged: any 
two nuclei among the TV spins significantly coupled to 
the electron are coupled to each other with comparable 
strength. The role played by this interaction in electron 
spin dephasing will be discussed in Section 6. Let us also 
mention, that recent experiments |50] are suggesting that 
this interaction is also affecting the dynamics of nuclear 
polarization (nuclear spin diffusion) at timescales much 
longer than the electron coherence time. 



3. Experiments in which the electron spin coherence 
is measured 

The spin dephasing affects all the experiments involv- 
ing coherent manipulation of the electron spin. We will 
outline now the basic experimental procedures used to 
gain information on decoherence of the electron spins. 
This will, however, be a theorist's description, in which 
we will assume that the electron spin is initialized (say 
in "up" direction), then well defined perfect rotations of 
this spin are driven by external stimuli at prescribed mo- 
ments of time, and finally the measurement (along a given 
axis) is done. Performing each of these steps for spins in 
quantum dots was in fact a major achievement of the last 
10 years, involving feats of experimental physics in the 
lab, and often requiring new theoretical ideas. Since we 
cannot give justice to these topics here, we refer the in- 
terested reader to review papers: about electrically con- 
trolled spins in gated QDs one can learn more in [4] (for 
the most recent experimental advances see e.g. [51 [TP] - ), 
while the optical control of spin rotations is described in 

[33137]. 
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Free Induction Decay 

Conceptually the simplest experiment is the Free In- 
duction Decay (FID), or simply the free evolution of the 
spin. We assume that the spin is initially oriented in, say, 
x direction (with B field, when nonzero, directed along 
the z axis). The spin is allowed to freely evolve for time 
t. At the final time t a measurement, say the projec- 
tion of the spin on the x axis, is made. If the spin was 
truly free (i.e. not interacting with any bath), its evolu- 
tion would simply be a precession, and the FID signal 
recorded as a function of the evolution time would be 
(&.(*)) = A cos nt. 

Such a simple oscillation is not observed for electron 
spins interacting with the nuclear bath. The reason is 
the spread of the Overhauser fields Vo in the ensemble of 
nuclear states, i.e. the inhomogeneous broadening. The 
currently made FID experiments are either spatial [51] 
or temporal ensemble [3J [7J [551 [S3J measurements. In the 
first case many spins (each interacting with a different 
nuclear bath) are measured simultaneously, while in the 
second case one measures repeatedly the same spin, but 
the data acquisition time is long enough for the nuclear 
bath to appreciably change its state. Assuming ergodic 
dynamics of the nuclei the two cases are equivalent, and 
they correspond to averaging of the electron spin evo- 
lution over the thermal ensemble of the nuclear states. 
This averaging will be done in Section 4. Here let us 
only say that it is enough to assume the nuclei static in 
order to obtain a very fast decay of the FID signal on the 
timescale of T| ~ y/N /A, which is about 10 ns in typical 
III-V dots. 

The T| time is an ensemble quantity. In a quantum 
computer one will deal with single qubits, and the single- 
spin decoherence time, T 2 , will be more relevant. Such a 
decay time would be observed if one could acquire enough 
signal from a single spin on a timescale shorter than the 
time on which the Overhauser field fluctuates apprecia- 
bly (this timescale is much longer than microseconds in 
GaAs dots pj [54]), or if one measured an ensemble of 
spins having the same Overhauser shift (as in |51j). In 
such cases one would be dealing with a so-called narrowed 
state of the nuclear ensemble. The free evolution experi- 
ment with such a nuclear ensemble is called the narrowed 
FID (NFID). 

The theory behind various ideas for achieving such a 
narrowed state is beyond the scope of this review (see 
e.g. [UJ and references therein). As for the experiments, 
the NFID measurement was done on an optically driven 
ensemble of InGaAs QDs [5T], in which the decay time 
of T 2 3 /is was seen. Recently a substantial narrow- 
ing of the nuclear state was also obtained in electrically 
controlled GaAs dots [55], leading to a tenfold increase 
(compared to T 2 ) of the decay time. A progress in single- 
shot measurement og GaAs spin qubits was also made [7J . 



Rabi Oscillations 

A two level system with level splitting ft, when exposed 
to a field of amplitude R coupling its levels and oscillating 
with frequency v, i.e. described by the Hamiltonian 

H = flS z + RS X cos vt , (8) 

will exhibit Rabi oscillations of the occupancy of its lev- 
els (see e.g. [55] )• These oscillations have the frequency 
il R = (R 2 + A 2 ) 1 / 2 , where A = ft - v is the detun- 
ing, and their amplitude is R 2 /Q 2 R . When the cou- 
pling between the system and its environment can be 
treated in Markovian approximation (i.e. when we can 
use Bloch equations), these oscillations are damped by 
exp(— t/T 2 ) factor. While this result is inapplicable to 
the case of the electron spin interacting with the nu- 
clear bath, the damping of the oscillations is of course 
expected. The Rabi oscillations were observed in electri- 
cally driven singlet-triplet GaAs qubit [3J , and in a single 
spin qubit [5J [351 HZ] ■ in [3H] the oscillations were visible 
for times up to a microsecond. A model explaining their 
decay will be briefly discussed in Section 4. 

Spin Echo 

The inhomogenous broadening, which obscures the in- 
teresting quantum dynamics of the single spin decoher- 
ence in the FID experiment, can be removed by a Spin 
Echo (SE) sequence [35] , in which the electron spin is ro- 
tated by angle 7r around one of the in-plane axes at the 
midpoint of its evolution. Such a protocol can be denoted 
as t/2 — 7r — t/2: two free evolution periods with a fast 
external pulse in the middle and readout at the final time 
t. It is easy to see that such a procedure will remove the 
static spread of the precession frequencies, since the evo- 
lution of every spin (assumed free, only having a random 
precession frequency) before the tt pulse will be undone 
after the pulse. This refocusing of the spins does not 
work perfectly when the bath is dynamic, so the ampli- 
tude of the SE signal is still decaying in time. In fact, if 
the Bloch equations were applicable, the SE would decay 
as exp(— t/T 2 ), with T 2 being the single spin dephasing 
time (the same would then be seen in NFID). In the more 
complicated case considered here, the decay of SE is non- 
exponential, and the characteristic time, Tse, needs not 
be related to the decay time of NFID. 

SE experiments were done in gated double GaAs dots, 
both for the singlet-triplet qubit [31 [HI HO] and for the sin- 
gle spin in one of the dots [5j . In the earlier experiments 
[3J [5] the magnetic field was B < 0.1 T, and the decay 
time was Tse < 1 /^s. This decay timescale was later 
explained by the theory [HI [20], which also predicted 
a characteristic oscillatory behavior of the SE signal at 
B > 0.1 T. This prediciton has been recently confirmed 
!9i], as discussed in more detail in Section 6. 
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Let us also mention that recently the SE experiment 
was performed with optically induced tt rotations on an 
ensemble of self-assembled quantum dots [llj . and on a 
single electron spin bound to a donor |581 , or confined in 
an InAs dot 1721. 



Dynamical Decoupling 

Application of a single 7r pulse during the qubit evolu- 
tion removes the effects of inhomogenous broadening, and 
furthermore it affects the dynamics of the whole system 
(qubit+nuclei). The application of multiple 7r pulses is 
well known in NMR |46j . where these pulses were shown 
to further suppress spin dephasing and relaxation. In 
the context of quantum computation these ideas were 
furhter developed under the name of dynamical decou- 
pling (DD), see e.g. [59 . For example, the Carr-Purcell- 
Meiboom-Gill (CPMG) sequence, which can be written 
asr/2 — 7T — t — 7T — — 7T — t/2 (with n pulses and the 
total sequence time t — nr, equivalent to SE for n — 1), 
was predicted theoretically to extend the coherence time 
of an electron spin in GaAs at high B fields [60] , and this 
prediction has been verified experimentally [9j. Other 
DD sequences were also considered theoretically for the 
case of spin decoherence due to hf interactions [2TJ1 IBD - 
|6"5] . Experimentally, in GaAs singlet-triplet qubits both 
CPMG [DJ EI] and other sequences with unequal spacing 
of pulses were employed [ID] . 



4. Quasi-static bath approximation 

As we mentioned before, the spin decay seen in the FID 
experiment can be calculated assuming completely static 
nuclei. This is equivalent to treating the '^2 i AiJi oper- 
ator as a classical field Bm- We neglect also the nuclear 
Zeeman energy and the dipolar interactions, and use the 
QSBA Hamiltonian Hqsba = ^S z + B n ■ S. For the un- 
polarized nuclear ensemble, in the limit of large number 
of nuclei N, one can use the Central Limit Theorem and 
derive the distribution of the Bm fields (see e.g. [2"Dll30| ): 



P(B 



1 



N ) 



(2^)3/2(73 



exp 



B N 

2a 2 



(9) 



with 



<7 2 = ^^ J Q (J« + 1)^A 2 = lj2 J ^ Ja + 1 ^ na 



A 2 
N 



where J a is the length of the nuclear spin of species a 
(e.g. J — 3/2 for all the nuclei in GaAs, and J = 9/2 for 
both isotopes of In). 

The calculation of the FID signal is especially easy 
at high B fields, when we can disregard the influence 
of transverse components of Bn and approximate the 



Hamiltonian by (f2 + B Z )S Z . Assuming that the spin is 
initialized along the x direction we have 



(Sx(«)) 



P{B N ) cos ((fi + B z N )t)dB N , 
cos(m)e- (t/T ^ 2 , (11) 



where T 2 * 



\[2/ a ^ \/N /A (e.g. for a bath with a single 
nuclear species and J =1/2 we have T| =^/8N /A). For 
typical III-V dots with N rj 10 6 nuclei we have T 2 * rj 10 
ns. The results for FID decay obtained within QSBA 
were confirmed by experiments 52. 53, 66, 67 . 

The calculation is slightly less trivial in the case of 
B = 0, where has to average the expression for the time- 
dependence of the spin in a classical static field 

S(t) = (So-n)n+(S — (So-n)n) cos B N t+S xn sin. B N t , 

^ (12) 

with the initial spin So and the direction of the effective 
field n = B N /B N , obtaining [55]I5D] 



{&(«)) 
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A plot of this result for GaAs with N = 10 6 nuclei is 
shown in Fig. [TJ Such a FID signal, which saturates 
at (S x (t)) = 1/6 for t > T 2 * was observed in 67. The 
maintaining of this saturation at long times is an artifact 
of QSBA, which has to break down eventually. Various 
theoretical approaches were used to predict the long-time 
1/ In t decay of the FID signal in this case [D5H7T] . 

While QSBA is surely valid on a timescale of T| , it is 
not known what is the precise limit on the timescale on 
which it is quantitatively valid (for a given experiment). 
One experiment, which sheds some light on this question, 
is the measurement of the decay of Rabi oscillations [3S] . 

In the rotating frame (in which the wavefunctions are 
transformed by the unitary operator exp(ii>S z t)), and 
within the Rotating Wave Approximation (in which the 
strongly oscillating terms are dropped from the Hamilto- 
nian), we have the effective QSBA Hamiltonian 



Hr,qsba = (A + B Z N )S Z +RS X 



(14) 



Assume that the spin is initialized in the "up" state at 
t = 0, and that the detuning A = 0. Furthermore, for 
simplicity let us assume R ^> B Z N (this does not affect 
the qualitative features of the result). The probability of 
finding the spin in the "down" state is given by [5pJ [3TJ 
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FIG. 1: Free Induction Decay (FID) of the S x component of 
the electron spin in a GaAs dot with N = 10 6 nuclei. The 
calculations are for B = and at B much larger than the 
rms of the Overhauser field, specifically B — 20 mT. These 
calculations are done using Eqs. (13 1 and (111, respectively. 



At longer times the B = signal is predicted to decay as 1/ lnt 
[68H71] . The FID decay signals such as shown in this Figure 
have been measured in many experiments [31 1521 1531 167] . 



The asymptotic ~ 1/yt decay predicted by the above 
formula was seen in GaAs spin qubit for times up to 
a microsecond [35], suggesting that QSBA is still valid 
at this timescale. The long-time 7r/4 phase shift of the 
oscillation was also visible in the data from [35]. The 
power-law decay of the Rabi oscillation was also seen in 
a Nitrogen Vacancy center spin qubit in diamond [72 . 

A calculation very similar to the one shown above can 
be performed for the case of the narrowed FID decay. 
Assuming for simplicity the narrowing condition B Z N = 
0, and assuming Bn <C fl (which corresponds to f2 3> 
a sa A/vN), by averaging Eq. (12 1 over the Gaussian 



distribution of Bfj and B V N we obtain an expression very 
similar to the previous one: 



(S*(t)) 



cos fit — rft sin fit cos (fit + arctan rjt) 



2(1 + V 2 t 2 ) 



(17) 

where we have defined r\ = o 2 /fl. Note that the charac- 
teristic decoherence time Tnfid, defined as the time in 
which the signal drops by half, is Tnfid ~ NO,/ A 2 , and 
it increases with increasing B field and dot size. 

In contrast to the Rabi oscillations discussed above, 
the envelope of the NFID signal decays asymptotically 
as 1/rjb, see Fig. [2] However, in the frame rotating with 
the fl frequency, the decay is (S x (t)) ~ I /rft 2 due to the 
long-time phase shift of 7r/2, which corresponds to the 
spin asymptotically rotating from the x axis to the y axis. 
This is an example of one of subtle effects associated with 
the bath-induced frequency shifts of the spin precession 
(see e.g. Ref. [5T| for discussion of other such effects in a 
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FIG. 2: The envelope of the Free Induction Decay for a nar- 
rowed nuclear state (NFID) in a GaAs dot with N = 10 6 
nuclei. The calculation is done using the QSBA (Eq. (17l). 



The results are shown for B = 30 mT and B = 100 mT. At 
higher B fields the decay becomes slower. At fields at which 
the decay timescale becomes longer than N/A~ 10 /xs, the 
QSBA becomes inapplicable, and the theories (both the ef- 
fective Hamiltonian based theory discussed in Section 6 and 
the Generalized Master Equation approach from 18, 21]) pre- 
dict a mostly exponential decay. 



different theoretical framework). 

In Section 6 we will show that Eq. (Jl7| can be derived 
as a short-time (t <C N/A) limit of a theory which in- 
cludes the nuclear dynamics. Since N/A is of the order 
of 10 fis in GaAs dots, this suggests that the QSBA is 
applicable at such defined short times. This agrees with 
the observation of a good agreement between the QSBA 
and the Rabi experiment for times up to a microsecond 



5. High magnetic field limit - dephasing due to 
dipolar interactions 

In large enough B fields one can completely neglect the 
flip-flop terms between the electron and the nuclei. The 
quantitative answer to the question how large these B 
fields need to be will be provided in Section 6. Now let 
us focus on the consequences of such an approximation. 
With Vg term neglected, we have the following Hamilto- 



H « Hz + i?dip + VqS z 



(18) 



The -ffdip is the only non-trivial term leading to nuclear 
dynamics. The crucial feature of this Hamiltonian is that 
it is of the pure dephasing form, i.e. the only electron spin 
operator present in it is S z . The obvious consequence is 
the fact that S z is then a constant of motion. A slightly 
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less obvious fact is the structural simplification of the 
theory of decoherence which takes place in this case. 

The physical picture of transverse spin dephasing is the 
following. The dipolar interaction leads to slow fluctua- 
tions of the Overhauser field: the nearby spins are flip- 
flopping because of the J^JJ term in Eq. M, and due to 
the spatial inhomogeneity of the hf couplings (i.e. Aj ^ Aj 
for the two spins considered) this corresponds to a change 
of the expectation value of Vq by Ay = Aj — Aj . This 
leads to fluctuations of the Overhauser field, and thus the 
fluctuations of the electron spin precession frequency. It 
should be noted that the hf interaction plays a passive 
role in this process: the nuclear fluctuations are due to 
the dipolar interaction, and the Overhauser operator sim- 
ply transmits these fluctuations to the electron spin. 

This process is known as spectral diffusion (i.e. the 
diffusion of the precession frequency of the spin). For 
many years it was treated using phenomenological clas- 
sical methods, i.e. by replacing the operator Vo by a 
classical stochastic process £oM ( see HU an( A references 
therein). The average transverse spin is then given by 



(S x (t))=Bc 



(19) 



where (...){ denotes the average over the realizations of 
the stochastic process £o(*)> and f(t;t') is a filter func- 
tion [2H1 [73] which parametrizes the sequence of tt pulses 
applied to the qubit (for FID f(t;t') — 1, and the filters 
for SE and 2-pulse CPMG are shown in Fig. |3| . Such an 
average can be easily performed when the process £0 (*) is 
assumed to be Gaussian (see e.g. [2H1I73]). A phenomeno- 
logical derivation of the properties of the £0 process was 
given in [25]. Let us also note that in the situa- 
tion when the decoherence occurs on times much longer 
than the typical interaction times between the spins in 
the bath (e.g. for dipolar coupling between the central 
spin and the bath spins), the assumption that £oW is 
an Ornstein-Uhlenbeck process was used to calculate the 
Rabi oscillations signal [74 , and it has been shown to 
describe very accurately experiments on dynamical de- 
coupling for Nitrogen Vacancy spin in diamond |75j . 

A fully quantum-mechanical solution to the spectral 
diffusion problem was given during the last 5 years [2"3T - 
[25] ■ Below we sketch this solution. 

Since the diagonal elements of the density matrix of the 
electron spin are constant in the pure dephasing case, we 
only need to consider the dynamics of the off-diagonal 
element p+_(t). The spin is interacting with the bath, 
so its density matrix is in fact a reduced one, obtained 
by tracing out the bath degrees of freedom from the total 
density matrix of the whole system (spin+ nuclei): 



! Tr B p(t) 



(20) 



where |±) are the eigenstates of the S z operator, and 
Ttb denotes tracing over the nuclear states. We as- 
sume a factorizable initial density matrix given by p(0) — 



PS (0)/3b (0) . For simplicity we will also assume that 
p± (0) = 1, and we will work with the decoherence function 
W(t)^p s ± (t)/P s ± (0). 

The operator U (t) describes the evolution of the whole 
system under the action of a series of ideal (<5-shaped) 
tt pulses applied to the electron spin. For pulses corre- 
sponding to rotations by angle tt about the x axis we 
have 



LT(t) = (-*)" e" 



iHr n 



o- x e 



-iHr n 



■ <J x e 



-iHri 



(21) 



with n being the number of applied pulses, Tj being time 
delays between the pulses, and the total evolution time 
t = Y^iii T i (SE sequence corresponds to n = 1 and t\ = 

T2 = */2). 

Plugging p(t) = U(t)p(0)U^(t) into Eq. (20) we arrive 



at the expression for the FID decoherence function 



WfidW 



e iH-t e -iH+t 



with 



H± — Hz + -ffdip ± - ^ MJt 



(22) 



(23) 



and (...) denoting the trace over the nuclei. Similarly, for 
SE we have 



Wsb(*) = U lH + T e iH - T e- lH + T e- lH - T 



(24) 



where r = t/2. Analogous expressions are obtained for 
DD sequences with more pulses. 



Equations ( 22 ) and ( 24 ) can be interpreted in the fol- 



lowing way: we have to calculate the evolution of the nu- 
clei under a time-dependent Hamiltonian (in which the 
sign of the Overhauser term is changing at the times at 
which the pulses are applied), first forward in time, and 
then backwards in time, but with a different Hamilto- 
nian. The presence of the electron spin is only causing 
H + ^H_. 

By going over to the interaction picture with respect 
to H^ (slightly complicated by the presence of tt pulses, 
see [20j). and after introducing the notion (well-known 
in the theory of nonequilibrium quantum dynamics, see 
e.g. [75]) of ordering of operators along a closed time- loop 
contour C, shown in Fig. |3j the expression for W(t) can 
be written in a very compact form 



W{t) = ( T c exp 



* J "Hdi P (T c )dr c 



(25) 



where 7c is the contour-ordering operator, and r c = (r, c), 
with t being the time variable and c= ± being the con- 
tour branch label (see Fig. [3]). The operators within the 
'Hdip('?c) interaction are given by 



J H T c) = ^ exp 



±i^r±i, j f{t;t')^dt' 



(26) 



c = + 



performing the trace. The result is 



T = 



b) 



-1... 



c) 



c = — 

f(t; t) for SE 



r = t 



T 



f{t; t) for 2-pulse CPMG 



r 



FIG. 3: a) The closed loop contour along which the operators 



in Eq. ( 25 1 are ordered, b) The plot of the time-domain filter 
function f(t; r) for the Spin Echo sequence, c) The same for 2- 
pulse CPMG sequence. The Figure is adapted from Ref. . 



For more details see |20j . The most important point is 
that in this Equation we have transformed our task into a 
calculation of an average of a generalized exponent (note 



also the similarity to the expression (19l). This can be 
done for any Hamiltonian of the pure dephasing type (see 
e.g. [77] for an application of this formalism to decoher- 
ence of a superconducting qubit). Such a task is very 
common in many areas of theoretical physics - a calcu- 
lation of a partition function of a given system being 
one example. Out of many methods devised for dealing 
with averaging generalized exponents, for our purpose 
the most useful is the linked cluster theorem (or the cu- 
mulant expansion) [751 EH] • According to this theorem 
W(t) can be written as an exponent of certain bath av- 
erages of products of Hdip operators: 



(27) 



where the v c terms (with the label c denoting sets of nu- 
clear indices) are such that each of them cannot be writ- 
ten as a product of averages involving subsets of nuclei 
from the set c. 

It turns out that in order to reliably calculate the NFID 
or SE signal due to the dipolarly-induced dynamics in 
III-V QDs it is enough to include only the second order 
i>2 terms (i.e. the terms involving the dynamics of pairs 
of nuclei) in the linked cluster expansion of Eq. ( |25| . 
This result has been obtained with a few methods |23F 
121)1 150] . The inclusion of higher order terms was shown 
to be necessary only when considering certain dynamical 
decoupling sequences [55] [53] . 

The formula for vi in the SE case is obtained by ex- 



v 2 



E 

i>j 



326^ 
A 2 



Aijt 



(28) 



where we have summed over the pairs of nuclei of the 
same isotope (the difference of Zeeman energies very 
strongly suppresses the dipolar interactions between het- 
eronuclear pairs at high B fields). This result can also 
be obtained by calculating the spin echo signal for only 
two nuclei coupled to the central spin, expanding the re- 
sult to the lowest order in bij, and summing over all the 
nuclear pairs. 

The above discussion was an attempt to outline the 
basic structure of the quantum theory of spectral diffu- 
sion. The reader interested in technical details is referred 
to the original papers [2"3Tf2"5] (among which [55] is prob- 
ably the most accessible for a newcomer to the subject, 
and the most detailed discussion of the diagrammatic 
linked cluster theorem for dipolar interactions is given in 
[80] ) . Here let us only state the results of calculations for 
GaAs dots. In these we have a large number N of rele- 



vant nuclei, and while each term in the sum in Eq. (28) 



panding Eq. ( 25 1 to the second order, using Eq. ( 26 ) , and 



is very small, their sum is large enough to cause decoher- 
ence on a timescale of ~ 10 /is (that is why it is enough 
to consider the pairs only). On this timescale we can 
use an approximation sinAjjt/8 ~ Ajji/8 in Eq. (28 1. 
This, together with Eq. ( |27| ), leads to prediction of the 
exp(— (i/TsE) 4 ) spin echo decay for GaAs. The charac- 
teristic time Tse = 10 — 100 /is depending on the shape 
and the size of the dot [24] [25] [27]. The NFID decay 
time was calculated to be of the same order of magni- 
tude [25] ■ A significant enhancement of the coherence 
time was predicted under dynamical decoupling [50lI63| . 
with the decay time prolonged from w 30 /is for SE to 
w 300 lis under a 6-pulse CPMG sequence in a GaAs QD 

As mentioned in the Introduction, these theories of SE 
decay due to spectral diffusion were successfully tested in 
measurements on spins of electrons bound to donors in 
Si. Very recently the SE decay was measured in a GaAs- 
based qubit for a wide range of magnetic fields [9J. The 
decay at B^>0.1 T was seen to be ^-independent, and to 
be well-fit by exp(— (£/Tse) 4 ) dependence with Tg E ~37 
/is. These observations are in a very good agreement with 
the theoretical predictions discussed above. 



6. Hyperfine-induced dynamics in moderate and low 
magnetic fields 

While the dipolarly-induced dynamics of the nuclei 
practically does not depend on the magnitude of the mag- 
netic field, the role of the processes involving the Vg part 
of the hf interaction is growing with decreasing B. This 
is most easily seen by looking at the f2 dependence of the 
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hf-mediated interaction derived in Section 2. At certain 
value of f2 it is expected that the electron spin decoher- 
ence will be dominated by purely hf-induced dynamics. 

Initial attempts [ST] to construct the decoherence the- 
ory employing the full hf Hamiltonian (and going be- 
yond the QSBA) were successful only in special cases 
(e.g. of nearly fully polarized bath), and standard time- 
dependent perturbation theory was shown to be of very 
limited use in this problem. Significant progress was 
made later by using the Generalized Master Equation 
approach j21j [32] and the equations of motion method 
for spin correlation functions [33J 131] ■ These approaches 
were used only in case of NFID, and assuming rather 
high magnetic fields (A/£l < 1 corresponding to B > 3 T 
in GaAs). 

The use of an approximate effective Hamiltonian con- 
taining the hf-mediated interactions allows for construc- 
tion of a theory which is structurally similar to the suc- 
cessful theory of spectral diffusion. We simply replace 
the dipolar interactions by the hf-mediated interaction 



19, 20 



in Eq. (25), with the only formal complication being the 
fact that the hf-mediated interaction from Eq. |7| is S z - 
conditioned (see [20] for details). 

There is one crucial difference with respect to the spec- 
tral diffusion theory. There, one only had to retain terms 
associated with small clusters of nuclei (e.g. pairs in the 
case of SE) in the cluster expansion of W(t). This was 
related to the fact that the dipolar interaction quickly 
decays with the distance between the nuclei, and at the 
timescale of the coherence decay the multi-spin nuclear 
correlations do not have the time to build up. The situa- 
tion is very different for hf-mediated interactions, which 
are very long-ranged, and one should include clusters of 
all possible sizes in the calculation. The assumption that 
only pairs contribute is valid only at high B fields (when 
A/Q < 1), at which such Pair Correlation Approxima- 
tion (PCA) calculations were first performed for the hf- 
mediated interactions [T7l [25l [61] . 



Luckily, there is a natural solution to this problem. 
When all the N spins are comparably coupled to each 
other, 1/N becomes a small parameter which controls 
the magnitude of contributions of various diagrams in the 
linked cluster expansion. It is possible then to sum all 
the diagrams of the leading order in 1/N, the so-called 
ring diagrams, and to obtain closed solutions for deco- 
herence under any sequence of pulses [19] [20] . Again, the 
formal feature of dealing with an average of a general- 
ized exponent is crucial: the solution is analogous to the 
calculation by 1/z expansion of the partition function of 
the Ising model with long-range interactions |82j (where 
z is the number of spins appreciably coupled with each 
other). 

Within this Ring Diagram Theory (RDT) one obtains 



W(t) w exp 



\k=2 



-i) k 



R k {t) 



(29) 



where Rk are the expressions for ring diagram averages. 
For NFID we have (assuming a single nuclear isotope for 
simplicity) 



Rk — 



J(J+1) 



3 

Af 



dA x . 



2 Al sin A i2 t sin A 23 t 



dA k p{Ax)...p{A k ) 
sin Ak\t 



(2fi)* 



A 



12 



.4 



2:S 



Ifcl 



(30) 



where p(A) is the distribution of the hf couplings. At 
short times t -C N/A we get Rk{t) {nt) k , and the 
summation of the series in Eq. ( [29] ) gives us the pre- 
viously obtained QSBA solution from Eq. (17). On the 



other hand, at long times 1 1$> N / A we obtain an expo- 
nential decay of NFID signal, W(t) « exp(— t/T 2 ), with 
T 2 ~NSl 2 /A 3 for A/Q<1 (this was also obtained in [T7] 
and [IB])- At these high B fields this solution agrees with 
the one obtained by a very different theory using the full 
hf Hamiltonian [21 . Thus, using the effective Hamilto- 
nian and the RDT we can reproduce the QSBA at short 
times, and also the long-time exponential decay which 
becomes important at high B fields (the 1/t 2 decay at 
even longer times predicted in [3TJ [331 131] is however not 
recovered by the RDT). 

The earlier PCA theory of decoherence due to hf- 
mediated interactions [TT] [5S1 [37] closely corresponds to 
retaining only the R 2 term in the sum in Eq. ( 29 ) . This 



is correct at high B fields (Q^>A), at which the calcula- 
tions in these papers were performed. 

For SE decay the RDT gave very characteristic predic- 
tions [19, 20 . At very high B fields the hf-mediated inter- 
actions between nuclei of different species are suppressed 
due to the differences Zeeman energies ui a p =u a — up. If 
one assumes a complete suppression of such inter-species 
interactions, the result is the absence of any decay of the 
SE signal [25]. This can be quickly checked by plugging 
the Hamiltonian with only intra-species hf-mediated in- 
teraction into Eq. (24). Since such an interaction term 



commutes with the nuclear Zeeman term, the four expo- 
nents cancel each other. In order to obtain the SE decay 
one has to sum all the ring diagrams which contain only 
the inter-species interactions. The result is particularly 
simple in the short time regime (t<^N/A): 



w SB (t) = 



i 

1 + R(t) 



with the function R(t) given by 



R(t) = 



AA 2 A 2 

N 2 n 2 uj 2 R 



a a apn a np sin 



4 U a pt 



(31) 



(32) 
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FIG. 4: Spin echo coherence function Wse in GaAs QD with 
N = 10 5 and 10 6 nuclei. The dots are calcul ate d usin g an 
approximate analytical expression from Eqs. (311 and (32 1, 



while lines are the result of the numerical summation of the 
ring diagrams (see [20] for details). The SE decay due to 
spectral diffusion (see Section 5) occurs on timescale of TJj| > 
10 us for these dots, and the decay in a SE real dot is given by 
a product of the results from this figure and exp[— (i/Isi) 4 ] 
factor. Such a decay was observed in recent experiments [9]. 
The figure is adapted from Ref. [19] . 



Finally, let us touch upon the question of limits of 
applicability of the RDT. Based on the analysis of the 
higher-order terms in the effective Hamiltonian, in [20] 
it was conjectured that this theory is reliable for mag- 
netic fields much larger than the rms of the Overhauser 
field (i.e. for 3> A/VN , which for a typical GaAs dot 
translates into B much larger than a few mT), at least 
at short times t <C N/A (corresponding to t <C 10 /is in 
GaAs). This conjecture can be considered experimentally 
confirmed in the case of SE [5] . The accuracy of the RDT 
for the SE calculation at low B fields was also confirmed 
by a study in which the RDT results were compared with 
the exact numerical simulations for a system of = 20 
spins [83] . where it was shown that general predictions 
of the RDT (e.g. a qualitative difference between hetero- 
and homo- nuclear baths) hold even at tt&A/y/N. 

At much higher fields (when 3> A) the RDT agrees 
with other theories [T5J [2T] at much longer times. The 
fact that in t <C N/A limit the RDT corresponds to the 
QSBA in the case of NFID, together with the success 
of QSBA in description of Rabi oscillation experiments 
at these times, supports our statement that low-i? and 
short-time behavior is correctly captured by the RDT. 



where a Q = |j Q (J a + 1). 

The SE signals calculated using the RDT are shown in 
Fig.[4]for GaAs dots at different magnetic fields. At high 
B, the hf-mediated interactions are causing the SE signal 
to oscillate with frequencies related to differences of nu- 
clear Zeeman frequencies. Note that because are not 
commensurate, the signal is not strictly periodic, and its 
revivals are never complete. With decreasing B the am- 
plitude of these oscillations grows, and the characteristic 
period becomes longer. Below the B c field corresponding 
to the electron Zeeman splitting f2 c « y/rA/N (where r 
is the ratio of electronic and nuclear Zeeman energies, 
r Ri 10 3 ), the SE signal goes practically to zero on the 
timescale of Tse ~ Sy/rT^ , which in GaAs translates to 
Tse ~ 100T 2 *. The partial revival of the signal at long 
times is then suppressed by the decay due to spectral 
diffusion, and this Tse can be considered a time of low-i? 
irreversible decay of the Spin Echo. 

The decay of the SE signal at such a timescale at low 
B was observed in earlier experiments [3] [6j. Very re- 
cently, the SE was measured in a GaAs singlet-triplet 
qubit for a wide range of B fields [H] • The characteristic 
oscillations and their evolution with B field was clearly 
observed there, and a very good fit of the RDT calcu- 
lation to these results was shown. At B > 0.3 T the 
SE decay due to the hf-mediated interactions was sup- 
pressed, and the decoherence was occurring due to the 
spectral diffusion. 



7. Summary and conclusions 

The dephasing of an electron spin by the nuclear bath 
is a nontrivial theoretical problem because of the strong 
(relative to the Zeeman energies at low magnetic fields) 
hyperfine coupling between the electron and the nuclei, 
and because of the slow dynamics of the nuclei. In fact, 
many experiments (Free Induction Decay for example) 
can be explained by a theory assuming that the nuclei 
are static. After a preparation of a special "narrowed" 
state of nuclei (in which the inhomogeneous broadening is 
decreased), or in a Spin Echo experiment, the coherence 
decay occurs at longer timescales, and the dynamics of 
the nuclei has to be considered. At high magnetic fields 
(larger than a few hundreds of mT in typical GaAs dots) 
the nuclei are fluctuating due to the dipolar interactions 
between them. The quantum theory of decoherence due 
to this process (the so-called spectral diffusion) success- 
fully explains the results of Spin Echo experiments at 
high fields (2"3T - f2l)] . At lower fields (between ~ 10 and 
a few hundreds of mT in GaAs) the hyperfine-mediated 
interactions among the nuclei are the dominant source of 
the nuclear dynamics causing the electron spin dephas- 
ing. The theoretical approach to this problem had given 
predictions for low-field Spin Echo decay [TjJ] [3U] , which 
have been recently confirmed experimentally [5]. 
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